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Programmable Matter Final Report 
Liu Team (Harvard) 


Project overview 


We propose the development and initial implementation of programmable matter, a new concept 
in materials. Programmable matter is constructed of mesoscale (mm-cm) building blocks that undergo 
self-assembly when exposed to a set of instructions specifying a desired three-dimensional target 
structure. When the target structure is no longer needed, the user can induce its efficient disassembly 
into building blocks that are suitable for future reassembly. 


To realize this ambitious concept, we propose the creation of polyhedral mesoscale building blocks 
("smartblocks") that are densely functionalized with DNA oligonucleotides on their surfaces. Each face 
of each smartblock is linked to many copies of a unique oligonucleotide, such that no two smartblock 
faces display the same DNA sequence. To induce the programmable assembly of these smartblocks 
into target structures, we propose to use "megastaples," long DNA duplexes terminated in 
single-stranded regions that are complementary to the oligonucleotides attached to two smartblock 
faces. Each megastaple therefore mediates the association of one pair of smartblocks via one 
smartblock-megastaple-smartblock interface, and a set of megastaples in principle can program the 
self-assembly of any stable three-dimensional structure. (Figure 1) 
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Figure 1 


The successful creation of programmable matter will require the development and integration of a 
number of concepts— some unprecedented— in computer science, chemistry, and materials science. 
We previously developed a framework for the development of these concepts in the form of three Tasks 
comprising nine Milestones. The Milestones for Phase I and Phase II are as follows: 

Milestone la: Develop a computational algorithm that outputs sets of hundreds to thousands of 
orthogonal ~50-base DNA sequences suitable for linkage to mesoscale polyhedra and suitable for the 
termini of megastaples. 

Milestone lb: Develop a software interface that accepts standard representations of three-dimensional 
structures together with a set of smartblock sequences and outputs a list of megastaple sequences 
suitable for mediating the self-assembly of the structure. 

Milestone 2a: Prepare and characterize a variety of candidate surfaces covalently linked to DNA strands 
of varying composition and linkage chemistries. 
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Milestone 2b: Develop a format for megastaples, experimentally validate the principle of DNA base 
pairing-mediated mesoscale surface assembly, and measure the interaction strengths afforded by such 
an approach. 

Summary of Progress: 

We have already achieved all three of the Phase I milestones listed above. We have developed 
and implemented an optimized algorithm that outputs hundreds to thousands of DNA 50-mer 
sequences that satisfy the orthogonality requirements that will be essential for our vision of 
programmable matter. We have also developed an interface with Visual OMP, an advanced 
commercially available oligonucleotide modeling platform, that enables the rapid screening of these 
algorithmically generated sequences to eliminate high melting-temperature spurious interaction. The 
generation, modeling, and screening of DNA sequences suitable for our early experimental efforts 
represent crucial capabilities that are now in hand (Milestone la). 

We also implemented in the form of a separate software program the ability to input a standard 
three-dimensional representation of an object (in the form of an AutoCAD file) and rasterize that object 
into smartblocks of arbitrary shape. We have tested this ability using basic shapes such as spheres, as 
well as more sophisticated arbitrary shapes such as wrenches. We are initially focusing our efforts on 
two smartblock shapes: cubes and truncated octahedrons. We believe these two smartblocks 
represent an attractive balance between ease of preparation and ease of self-assembly. 

We have also developed and tested a variety of surface functionalization chemistries that enable 
the covalent attachment of DNA to glass and PDMS surfaces. We implemented the most promising of 
these chemistries and measured the functionalization density on glass and PDMS surfaces as ~10 12 and 
~10 13 DNA 50-mers per 0.25 cm 2 . Both functionalization densities are already well above our stated 
goal of at least 10 11 DNA 50-mers per 0.25 cm 2 of surface. 

In addition, we completed a portion of Milestone 2b during Phase I, by synthesizing and 
characterizing megastaples. We used a PCR-based method to synthesize the megastaples using 
primers attached to our handle sequences using an unnatural spacer. We then confirmed the 
structure of our megastaples through both gel electrophoresis, and a fluorescence binding assay. We 
also conducted extensive experiments on three different approaches to DNA hybridization-mediated 
surface-surface interaction. Most recently, we have embarked on a new approach to achieving 
DNA-mediated surface associations that uses extremely long, highly repetitive single-stranded DNAs 
linked to a variety of surfaces. This new approach in principle offers much more flexibility for 
productive DNA hybridization and potential interaction strength than our initial approach, and is 
currently being tested for its ability to mediate detectable sequence-specific surface interactions. 

We made substantial Phase II progress as well. Specifically, we used surface chemistries described 
below to generate DNA-linked glass and PDMS surfaces, and demonstrated surface-surface association 
that is dependent on DNA sequence, DNA structure, and DNA hybridization. Finally, we also applied 
these principles to the assembly of DNA-programmed plasmonic nanoclusters of gold spheres and gold 
shells functionalized the DNA oligonucleotides. 
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Milestone la: Sequence generation 


Overview 

J-Seqgen is software which generates a series of DNA sequences which are orthogonal to each 
other, i.e., they should not form complexes with each other above a certain threshold temperature. 

The input parameters are: strand length, number of strands, criterion length, and threshold 
temperature. The program also has the ability to preload sequences to be avoided during the main 
sequence generation. The program can output a human-readable text file and/or an OMP-compatible 
file containing the sequences. 

The program generates the sequence using sequence symmetry minimization (SSM), and then 
optionally checks the melting temperatures using the commercial software Oligonucleotide Modeling 
Platform Developer Edition (OMP) from DNA Software. 

Background 

Sequence symmetry minimization: 

SSM is a heuristic method that designs a sequence with no unintentional spurious interactions 
of a specified length. The essence of SSM is that, if a subsequence of a certain length appears in the 
sequence, its compliment should not appear anywhere in any of the sequences that it is not designed to 
complex with. The subsequence length used in imposing this requirement is called the criterion length 
and each subsequence of this length is a word. Sequences which fulfill SSM will have no spurious 
complimentary runs of at least the criterion length in any frame. For example, if the criterion length is 
7, there will be no possibility of forming spurious duplexes of 7 or more bases. 

Using a lower criterion length lowers the strength of undesired interactions; however, the size 
of the sequence set to be generated places a lower limit on the criterion length. For criterion length c 
there are only 4 C words available to be used. In our system, with m strands of length n plus their 
antistrands, the total number of words contained in a single strand is n-c+1, and the number of words 
contained in the whole system is 2m(n-c+l). Furthermore, the user may wish to exclude certain words 
such as those contained in the megastaple template. Since each word must be unique, c must be 
chosen so that 4 c -p > 2m(n-c+l), if p words are excluded. 

Oligonucleotide modeling platform: 

OMP is commercial software which calculates thermodynamic parameters such as the free energy 
change of hybridization (AG) and the melting temperature (T m ) for collections of strands. J-Seqgen 
interacts with OMP through the functionality provided in the Developer Edition of the software. 
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The algorithm 


Pseudocode: 

Preload excluded sequences into word hashtable. 

For each strand: 

Find an unused word to start the sequence. 

While unassigned bases remain: 

Check whether adding each of the four bases is allowed. 

If no base is allowed, back up one base. 

Otherwise pick a random base from those allowed, advance to next base. 

Once a candidate strand is generated, use OMP to check T m of all heterodimers with previously 
generated strands; reject and regenerate if any T m ’s are greater than user-defined threshold. 

Detailed explanation: 

J-Seqgen uses a hashtable to store whether each of the 4 C possible words has been previously used 
and thus is unavailable to be used again. The advantage of a hashtable is that it can be accessed in 
constant time, which is vital for the linear scaling of the algorithm. 

The program's user interface allows the user to preload words into the hashtable from a file, 
set the generation parameters (strand length, number of strands, criterion length, and melting 
temperature threshold), start the sequence generation, and output an OMP OEF file. 

The sequence generation itself function starts by generating a random word to use as the first 
word in the strand; if this word has already been used, another random word is tried. Once an 
acceptable word is found, the program starts adding new bases one at a time. All four words 
generated by adding each of the four nucleotides are checked in the word hashtable. If a word has 
already been used, then that base is rejected. If all bases have been rejected, the loop backs up to the 
previous base. Otherwise, the algorithm randomly chooses one of the acceptable bases, and moves to 
the next base. 

Once one complete strand sequence has been generated, the program then checks the 
candidate strand with OMP. It checks the pairwise melting temperature of the candidate strand with 
each of the previously generated strands, and considers the maximum of these. If this maximum 
melting temperature exceeds the user's threshold temperature, the algorithm rejects the strand, and 
generates a new one. 

Running time: 

The SSM algorithm runs in nearly linear time, since one extra word needs to be checked for each 
extra base to be generated. The algorithm deviates from linear time in that long sequences will back 
up more, but this is not observed to be significant unless the number of needed words is very close to 
the number of possible words for that criterion length (4 C - p). 

OMP is observed to run in linear time with respect to the number of pairwise interactions. This 
gives an overall complexity of 0(m 2 ), where m is the number of strands. The complexity with respect 
to the strand length n has not been thoroughly examined but it seems to be much worse. 
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Figure 2. Histograms of the melting temperatures of pairwise spurious interactions in three 
different sets of 400 50-mers and their complements. The range for each data point in the 
histogram is 1 C°. Top: Sequence set generated using SSM only. Middle: Sequence set 
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generated using SSM, and using a threshold temperature of 50°C. Bottom: Sequence set 
generated using SSM, and using a threshold temperature of 45°C. 

Overall, 1000 50-mers can be generated in <1 minute using only SSM. If OMP is used to check 
melting temperatures without rejecting any strands, the same generation will take =12 hours. Using a 
temperature threshold will increase the time needed since strand must be rejected and regenerated; a 
threshold temperature of 50°C will take =24 hours and a threshold temperature of 45°C will take nearly 
a week. 

Results 

Programming: 

Melting temperature data for sets of 400 50-mers and their complements was generated (Figure 
2). This comprises 800 strands total and 319,200 interactions. These were generated using sequence 
symmetry minimization (SSM) algorithm, both with no temperature threshold, and with checking by 
OMP using a threshold temperature of 50°C or 45°C. 

Empirical confirmation: 

Gel electrophoresis was performed (Figure 3) to demonstrate that duplexes only form between 
the intended sense and antisense handle strands generated by our software. Four sets of strands were 
arbitrarily chosen from our 1000-strand collection and incubated at 50°C in various combinations. 

Bands representing duplexes are visible in lanes where duplex formation is intended; no bands 
representing unintended interactions are observed in any lane in the gel. 
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Figure 3. Gel electrophoresis of sequences. A prime symbol indicates the "antisense" 
complement of a strand. The left half of the gel shows single strands and duplexes, followed by 
combinations of all sense, all antisense, and all strands. The right half of the gel shows strand 
999 plus each of the other strands; no interactions are observed except for the desired 
interaction with its complement. The experimental conditions were 1 pmol ea. strand in PBS 
buffer, incubated at 50°C for 1 hr, 2% agarose+EtBr, 100V, 25min. 


6 

























Milestone lb: Object rasterization 


Overview 

J-Honey is software which rasterizes an arbitrary shape into a set of polyhedra (smartblocks) 
which approximates it. It takes two inputs: the honeycomb lattice type, and a function which specifies 
whether a given set of coordinates is inside the desired shape. The output is a human-readable file 
containing the smartblock specifications, and/or an AutoCAD script representing the assembled 
smartblocks. 

The smartblock specifications consist of, for each polyhedron, a polyhedron ID and the 
polyhedron type; and for each face, a sequence ID as well as the polyhedron ID of the complementary 
smartblock. 


Background 


A honeycomb lattices is a three-dimensional tessellation of solid polyhedra which completely fills space. 
There are 28 ways to fill space with convex uniform polyhedra: 


cubic 

rectified cubic 
truncated cubic 
cantellated cubic 
bitruncated cubic 
cantitruncated cubic 
runcitruncated cubic 
omnitruncated cubic 
alternated cubic 
truncated alternated cubic 
runcinated alternated cubic 
cantitruncated alternated cubic 
quarter cubic 
gyrated alternated cubic 


gyroelongated alternated cubic 
elongated alternated cubic 
gyrated triangular prismatic 
gyroelongated triangular prismatic 
truncated square prismatic 
snub square prismatic 
hexagonal prismatic 
truncated hexagonal prismatic 
trihexagonal prismatic 
triangular prismatic 
rhombi-trihexagonal prismatic 
omnitruncated trihexagonal prismatic 
snub trihexagonal prismatic 
elongated triangular prismatic 


5 of these use only one type of 
polyhedron: 

cubic: cubes 

bitruncated cubic: 

truncated octahedra 

triangular prismatic: 

triangular prisms 

gyrated triangular prismatic: 

triangular prisms 

hexagonal prismatic: 

hexagonal prisms 


There are many other ways to fill space with a single type of non-uniform polyhedron. The program 
currently works with the cubic (cubes) and bitruncated cubic (truncated octahedra) honeycombs. It 
can be upgraded in future versions to work with more or custom honeycomb types. 
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The algorithm 


Pseudocode: 

Initialize: get honeycomb parameters from user 

Create first polyhedron 

For each face of current polyhedron 

Determine coordinates of adjacent polyhedron. 

Query whether coordinates are inside object. If no, discard new polyhedron. 

Query whether coordinates are already occupied. If yes, assign complimentary sequence ID's 
to corresponding faces. 

Otherwise, create new polyhedron and place at end of queue. 

Rerun for-loop on next shape in queue. 

Detailed explanation: 

The honeycomb type and a function describing the shape to be rasterized are set by the user. The 
JHoney object contains a queue of polyhedra stored in a linked list. The algorithm creates the first 
polyhedron at coordinates 0,0,0 and then considers each shape in the queue until it reaches the end. 

For each polyhedron, each face is considered individually in turn. 

For a each face of a polyhedron, the algorithm first determines the coordinates of the polyhedron 
adjacent to the current face. It then determines whether this point is inside the object. If it is 
outside, the face is assigned to be an exterior face and the function returns, otherwise the function 
continues. The algorithm then checks if there is already a polyhedron at those coordinates. If yes, 
the sequence ID's of both complementary faces are set; if no, a new polyhedron is created with those 
coordinates (with the sequence ID's remaining unassigned), and placed at the end of the queue of 
polyhedra. 

Once the algorithm is finished with one face, it runs subsequent faces on the same polyhedron, and 
eventually the program advances to the next shape in the queue, until the end of the queue is reached. 
Data can be output as either a human-readable text file or an AutoCAD script. 

JHoney4ARX is additional code which allows the J-Honey software to be accessed from within the 
AutoCAD evironment. It uses the ObjectARX for AutoCAD 2008 code library available from the 
Autodesk website. This code produces a file with a .arx extension. This file can be loaded into 
AutoCAD to activate the commands. JHoney4ARX allows the user to rasterize the object input into the 
AutoCAD program and output the human-readable or script file. 

Results 

Sample output is provided (Figure 4). Arbitrarily-shaped objects in AutoCAD can be rasterized into 
either cubes or truncated octahedra, along with the specifications for assigning DNA sequences to each 
of the smartblock faces. 
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Figure 4. This demonstration of the Milestone lb software shows an arbitrarily-shaped object 
rasterized into different honeycomb lattices. Top: the original objects, a sphere and a wrench. 
Middle: the objects as approximated by cubes. Bottom: the objects as approximated by 
truncated octahedra. This figure demonstrates both the success of our honeycombing 
algorithm as well as our ability to use an arbitrary, irregular shape from AutoCAD as an input. 
The honeycombing algorithm also generates a DNA sequence identifier for each face of each 
polyhedron, not shown. 


9 







Milestone 2a: Surface functionalization 


Overview 

Milestone 2a requires us to develop the surface chemistry needed to attach DNA oligonucleotides 
to a PDMS surface. We investigated a number of chemistries and characterized their suitability for 
DNA surface functionalization. 

Results 

PDMS oxidation: 


The first step of surface functionalization is oxidation of the PDMS surface. Untreated PDMS 
surface is very hydrophobic. After plasma treatment, the surface becomes much more hydrophilic due 
to the formation of polar hydroxyl groups on the surface. We have found a good set of conditions at 
which we can obtain a hydrophilic PDMS surface (see Figure 5). 
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Figure 5. Top: The chemical scheme for chemically functionalizing glass and PDMS surfaces. 

Bottom: A 3 pL water droplet on different PDMS surfaces. We can see that after a 50 mA, 10 min 
oxygen plasma treatment, the PDMS surface became rather hydrophilic and water easily spread 
out, which is the desired result. After silane deposition, the surface becomes hydrophobic 
again, as expected. 

Coupling reactions: 

We tested four different DNA coupling reactions (Figure 6). The purpose of testing these 
reactions is to select the best one for the large-scale fabrication. 
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Figure 6. The four DNA coupling chemistries we considered. 
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The coupling efficiency was measured using a fluorescence assay. A fluorescently-labeled DNA 
oligomer was hybridized to the DNA immobilized on the surface, and the fluorescence was measured 
(Figure 7). 
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Figure 7. Fluorescence assay to evaluate the efficiency of the four coupling reactions we 
considered. In each of these images, the left column is the experiment and the right column is 
a negative control. 


We selected the epoxide-amine reaction as the best candidate reaction based upon the coupling 
efficiency as well as its procedural simplicity and the availability of materials. Therefore, we decided to 
choose this reaction for all the following experiments. 


An experiment was designed and carried out to quantitate the surface DNA density (Figure 8). 
We directly deposited different amounts of Cy3-labeled DNA on a substrate (either glass or PDMS), in 
order to create a working curve relating fluorescence intensity to fluorophore density. We then used 
this working curve to quantitate the density of a DNA-modified surface, after hybridization with 
Cy3-DNA, based on the fluorescence intensities. 
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Figure 8. Covalently attached fluorescently-labeled DNA is compared to standards of known 
concentrations of free Cy3. Left: glass surface; right: PDMS surface. Plotting the DNA 
fluorescence on the standard curve yields DNA densities of 3.2xl0 12 and 3.3xl0 12 
molecules/0.25cm 2 respectively, surpassing our original goal. 

These results indicate that we have achieved the goal (10 11 molecules/0.25 cm 2 ) from our 
original proposal. Our DNA density after hybridization is greater than 10 12 molecules/0.25 cm 2 . 
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Milestone 2d: Megastaple Synthesis 


Overview 

Milestone 2d requires us to develop a format for megastaples, which are about several thousands 
base pairs long DNA duplexes with two 50-mer single stranded overhangs on the both termini. 

We have designed our megastaples as 2 kbp duplex DNA's with single-stranded overhangs 
complimentary to the smartblock handle oligonucleotides. This is for two reasons. First, the long 
duplex region will bridge micron-sized microscopic imperfections in the polyhedron surfaces (Figure 9). 
Second, it allows for all the assembly information to be encoded in the megastaples, so that different 
sets of megastaples will result in different controlled assemblies using the same smartblocks. 
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Figure 9. Schematic of the target Milestone 2b setup. Two surfaces are functionalized with 
oligonucleotide handles. Upon addition of megastaples containing overhangs complimentary to 
the surface handles, the two surfaces will adhere to each other. Addition of noncomplimentary 
megastaples would not result in adhesion, providing the selectivity vital for our proposed 
self-assembly method. Additionally, the length of the megastaples would allow them to bridge 
micron-scale imperfections in the surfaces. 


Experimental design 


We have used a polymerase chain reaction (PCR)-based method to synthesize megastaples. Our 
method is similar to standard PCR using a plasmid as a template, except that we use special primers 
which are connected to our handles via an organic spacer. The spacer blocks the polymerase from 
acting on the handle portion of the primer complex, which results in double-stranded megastaples with 
the handles remaining as single-stranded overhangs. This allows us to efficiently create the 
megastaple's duplex portion while retaining full control over the handle sequences. 

It is critical that the handle regions of the megastaples remain single-stranded, since these comprise 
the sticky ends that will interface with the smartblock handles. We can confirm that the 
single-stranded regions exist through their ability to bind fluorescently-labeled oligonucleotides which 
are the compliments of our handle sequence (Figure 10). If the megastaple has the proper structure, it 
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should become fluorescent once annealed to the handle compliments, which can be confirmed by gel 
electrophoresis. 


spacer 



_ Flu 

Cy3 


Figure 10. Scheme for synthesis and structure confirmation of megastaples. A polymerase 
chain reaction (PCR) is performed using special primers which are connected to our handles via 
an organic spacer, blocking the polymerase from acting on the handle region. The existence of 
the single-stranded handle regions can be experimentally determined by annealing the 
megastaples with the fluorescently-labeled handle compliments, and performing gel 
electrophoresis to determine megastaple fluorescence. 


Results 

We have verified the structure of the new megastaples by performing the fluorescence 
hybridization experiment explained above. The megastaple is clearly labeled by both handle 
compliments, meaning that sticky ends are available for smartblock hybridization. (Figure 11) 


o 

E 



Figure 11. Images of a gel electrophoresis analysis of the fluorescently-labeled megastaple 
structure in Figure 10. a: Cy3 fluorescence, b: Fluorescein fluorescence c; Overlay of Cy3 
and fluorescein fluorescence, d: Fluorescence after staining with ethidium bromide. The gel 
conditions are: 1% agarose, 100 V, 28 min. The lane labels over image c are valid for all images. 
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Based on a similar strategy, megastaples with length up to 8k base pairs were also successfully 
synthesized. The reason to develop longer megastaples is that surface defect is one of the major 
problems that inhibit DNA-mediated surface interaction. The 8k base pairs megastaples are expected to 
be able to tolerate surface defects up to ~2.5um, thus facilitating the surface interaction. Our PCR 
results showed the long megastaple synthesis was very successful (Figure 12). 
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without template (lambda DNA) 

Figure 12. Agarose gel electrophoresis analysis of the PCR products. The desired length of the 
megastaple is about 8k base pairs, as indicated in the gel. 
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Milestone 2b: DNA-Mediated Surface Interaction 


Overview 


The target of Milestone 2b is to experimentally validate the ability of DNA hybridization to mediate 
the association of two DNA-conjugated surfaces. This is the key step of the entire project. Using 
molecular recognition to direct the assembly of mesoscale objects is a new fundamental concept in 
basic science and hasn't been demonstrated yet. Our recent efforts to validate this concept include two 
different approaches: (1) Using megastaples to direct the association of the mesoscale objects; (2) Using 
extremely long repetitive single-stranded DNA to direct mesoscale surface association. 

(1) Megastaple-Directed Mesoscale Surface Association 

Experimental design 


We developed a fluorescence-transfer experiment to characterize the efficiency of 
megastaple-mediated surface interactions. In this experiment, fluorescent probe labeled megastaple is 
first hybridized to the surface-attached single stranded DNA. Then another surface with the 
complementary DNA single strands is brought into contact with the original surface. If the hybridization 
takes place, after separate the two surfaces, we should be able to see the fluorescence signals on the 
second surface (Figure 13). 
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Figure 13. Scheme for a fluorescence-transfer experiment to characterize the efficiency of 
megastaple-mediated surface interactions. Fluorescent megastaple is first hybridized to 
oligonucleotides immobilized on Surface 1. Surface 2 is then brought into contact with Surface 
1. After the two surfaces are separated, fluorescence should be observed on Surface 2 only if 
the megastaples were able to hybridize both surfaces simultaneously. 


Results 

From the experimental results, we can see that after surface separation, there is some fluorescence 
signal showed up on the second surface (Surface 2). The surface 2 is PDMS with 100pm protrudent 
square patterns. In the image, only the patterned area, which could effectively contact the Surface 1, 


17 
























has signal, suggesting that the fluorescence transfer is due to a hybridization event involving the two 
surfaces in close contact. 



Figure 14. Fluorescence image of Surface 2 after DNA hybridization and surface separation. 
The fact that fluorescence signal is present on Surface 2 indicates that megastaples were able to 
hybridize both surfaces simultaneously. Additionally, the pattern of fluorescence corresponds to 
the raised areas patterned on to Surface 1 (microscope image, inset), suggesting that the 
fluorescence transfer is due to a hybridization event involving the two surfaces in close contact. 


However, there are two major problems existing in the experiment. First, the reproducibility of this 
experiment is not very high (~30% so far), and sometimes we also got very strong background signals. A 
possible reason is that, we couldn't wash off all the free DNA strands after DNA attachment and 
hybridization, probably due to the hydrophobicity of the PDMS surface. Secondly, after the hybridization 
takes place between the two surfaces, we couldn't observe very strong interaction between the surfaces, 
which indicates the hybridization yield is still very low. In order to improve the reproducibility and also 
hybridization yield, we also tried a number of different approaches. 

Approach 1: DNA-mediated interaction between gold surfaces 

We tried to develop DNA-mediated gold surface interactions. Gold surface has some distinct 
properties other than glass and PDMS. Since flat gold surface can be easily made by metal evaporation 
and the gold-thiol interaction is very reliable and easy to carry out, it's quite feasible to introduce single 
stranded DNA onto the gold surface. Besides, free DNA single strands can be easily cleaned off from the 
gold surface, thus effectively lowering the background signals. Our experimental results showed that the 
fluorescent probe DNA can effectively hybridize to DNA that is attached onto the gold surface, while the 
negative control sample gave very weak fluorescent signal, indicating a much lower background (Figure 
15a), as we expected. 

We also performed the gold surface interaction experiment, as shown in Figure 15b. The thiolated 
ssDNA was first immobilized onto the gold surface. The surface was then passivated with 
mercaptohexanol. Subsequently, another complementary thiolated ssDNA was added and hybridized to 
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the surface bounded ssDNA. After a thorough washing step, another bare gold surface was brought into 
contact with the DNA modified gold surface and the two surfaces were incubated for overnight. After 
incubation, we observed some weak interaction between the two gold surfaces. And the negative 
control sample, which doesn't contain the thiol group, didn't show any interaction between the 
surfaces. This result was quite promising. However, the interaction was still quite weak. We 
hypothesized that the rigidity of the gold surface would inhibit the good conformal contact between the 
two surfaces, thus there was only a small portion of the surface bounded DNA can effectively participate 
the hybridization and the surface reactions, resulting in a much weaker surface interaction. Therefore, a 
softer and squeezable surface is expected to be a better choice. 

a) b) 
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■ weak interaction observed 

I (sassixsto)) 

• M group 

Figure 15. a) Fluorescence image of the gold surface after hybridization with the fluorescent megastaples. 
Left, the surface bounded ssDNA is complementary with the fluorescent megastaples; Right, the surface 
bounded ssDNA is not complementary with the fluorescent megastaples, b) The scheme of testing the gold 
surface interaction mediated by DNA hybridization. 
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Approach 2: DNA-mediated interaction between dendrimer coated surfaces 

Dendrimers are repeatedly branched molecules, which could contain lots of hydrophilic functional 
groups. Using dendrimer coated surface as substrate to perform DNA-mediated surface interaction 
would have several potential advantages: 1) the densely packed functional groups could react with DNA 
molecules, thus increasing the surface DNA coverage and surface hydrophilicity. 2) Dendrimer molecules 
are about several nanometers in size. If the surfaces are modified with these bulky molecules, the 
surfaces would become softer and may better contact with each other, which could, therefore, facilitate 
DNA hybridization. 

The protocol of immobilizing dendrimer molecules onto the surface is illustrated in Figure 16a. 

After immobilization of dendrimer molecules, single stranded amine DNA was added onto the surface, 
followed by the addition of the complementary fluorescent probe DNA. The hybridization assay results 
on both glass substrate and PDMS substrate are shown in Figure 16b. From the results, there are a 
number of things we learned: 

(1) ssDNA can be covalently linked onto the surfaces and further hybridize with complementary 
fluorescent DNA. 
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(2) Dendrimer modification didn't improve the DNA density quite much compared with the 
non-dendrimer surface. 

(3) The glass surface, after modified with dendrimer and DNA, became fairly hydrophilic. However, 
the PDMS surface, after modified with dendrimer and DNA, was still not hydrophilic. This is a 
potential problem for the surface interaction. 

(4) The PDMS surface became a little blurry after all steps of modifications, probably due to the 
solvent-induced swelling and acid etching. Therefore, this protocol might not be a good way for 
the PDMS surface modification. 

a) 

OH OH OH OH APTS NH2NH2NH2 NH2 



(64 -COOH group* 

-» 

EDC, s-NHS 


Glutaric 

anhydride 




modified with dendrimer without dendrimer modified with dendrimer 


Figure 16. a) The procedure of immobilization of dendrimer molecules onto the surfaces, b) Fluorescence 
assay of the DNA hybridization on the dendrimer coated glass surface (left) and PDMS surface (right). 


Approach 3: DNA-mediated interaction between polymer-coated surfaces 

We also tried to immobilize DNA single strands onto the polymer coated glass slide. The polymer 
we chose has hydrophilic backbones. It also contains NHS ester that could react with amine DNA. The 
thickness of the polymer layer is about lOnm, which provides a slightly soft, squeezable surface. We 
used this polymer coated surface as substrate (Surface 1) and performed the fluorescence transfer 
experiment again. 

In order to improve the reliability of the fluorescence-transfer experiment, we have redesigned the 
experiment to use megastaples with two fluorophores directly attached to them. Compared with our 
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previous version of fluorecent labeled megastaple, the new one has sevearal additional advantage: 1) 
The fluorophore is covalently attached to the megastaple, meaning we do not have to worry about an 
extra intermolecular interaction; 2) One of the handles is no longer occluded for half its length. 3) Since 
there are fluorophores on both sides, there should be no difference in surface quenching on the top and 
bottom surfaces. 

The experiment design is shown in Figure 17b. From the results, we observed fairly strong 
fluorescence signal on the positive PDMS surfaces after the transfer experiment. Meanwhile, the 
negative PDMS samples showed much lower fluorescence signal. The polymer coated glass slide also 
showed reasonable pattern after the transfer. The fluorescence signal in the contact area became much 
weaker. Moreover, we also observed some interaction between the PDMS and glass surfaces. All these 
phenomena suggested the DNA hybridization could mediate mesoscale surface interaction. 

However, there are still a number of disadvantage of the polymer coated surface. Currently, the 
polymer coated surface is only available on the glass slides. It would be desired to prepare such a 
polymer coating on some other substrates. In addition, for the negative samples, we also observed some 
interaction between the two surfaces. That's probably because the polymer layer is relatively soft and 
squeezable. It could change surface morphology upon pressure during the DNA hybridization. This 
morphology change could also introduce some non-specific interaction between the surfaces. 

Therefore, the difference between the positive sample and the negative sample is not quite clear. We 
still need to further enhance the DNA hybridization induced specific surface interaction. 


21 



a) 


Surface 2 


Surface 2 
(PDMS> 


*\ 


A' 

+ 


Surface 1 
(polymer coated 
glass slide) 


Surface 2 




Surface 1 


B 


separate 



Surface 1 
(polymer coated 
glass slide) 



positive 

Surface 2 has 
DMA strand B 


v j 

negative 1 

Surface 2 has 
DNA strand B' 


negative 2 

Surface 2 has 
DNA strand A 


surface 2 


Figure 17. a) Scheme of a new version of fluorescence-transfer experiment. Note that in this version, a 
symmetric fluorescence labeled megastaple was employed, b) Fluorescence transfer experiment results. 


Research In Progress Towards Milestone 2b (and the other Phase II milestones) 


(1) Extremely long repetitive ssDNA-mediated surface association 


Although our initial efforts on the megastaple-mediated surface interaction showed promising 
results, it still remains to be confirmed thatDNA hybridization can induce strong mesoscale surface 
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interaction. In order to achieve more convincing results, we developed a new approach. Instead of 
using 50 base pairs of DNA hybridization to induce the surface interaction, we are now using 
single-stranded DNAs thousands of base pairs in length with highly repetitive sequences. If we can 
modify the two surfaces with complementary single stranded DNA with length of several thousand 
bases, the potential interaction between the two surfaces could be even higher than the theoretical 
Newton-scale interactions of the original design. As a result, realizing even a small fraction of the 
theoretically possible interactions in this new design may enable our basic vision of programmable 
matter to be realized. 


Experimental design 

We proposed to use rolling circle amplification (RCA) to prepare super long repetitive DNA single 
strand. The mechanism of RCA reaction is shown in Figure 18. It utilizes a small circular single stranded 
DNA as a template. At the presence of primer, polymerase, and dNTP, the primer can elongate along 
the circular template for multiple cycles, generating a super long (several thousand bases or even longer) 
repetitive single stranded DNA product. Another advantage of the repetitive sequences is that the two 
complementary DNA single strands can hybridize at any position. Consequently, using long repetitive 
ssDNA to mediate surface interaction could potentially tolerate much larger surface defects. 



polymerase 

dNTP 



Figure 18. Mechanism of the rolling circle amplification reaction. 


We proposed two approaches to carry out the experiments: 

1. Approach A: 

(1) Using RCA to generate two long repetitive single stranded DNA. The two primers both have a 5' 
amino group. The two circular templates are complementary to each other. 

(2) Immobilize the obtained single-stranded DNA onto two surfaces 

(3) Test and quantify the surface interactions. 
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2. Approach B: 

(1) Immobile the two aminated DNA primers onto two difference surfaces 

(2) RCA on surface to generate surface bounded long repetitive single-stranded DNA. 

(3) Test and quantify the surface interactions. 


Initial Results 


1. ligation of the circular template: 

In order to perform RCA, we need a circular DNA as template. Therefore, the first task is to 
synthesize circular ssDNA via ligation. The ligation result is shown below: 
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Figure 19. a) TBE-Urea polyacrylamide gel electrophoresis (PAGE) analysis of the ligated circular DNA 
templates, b) LC-MS characterization of the purified circular monomer DNA templates. 


The ligation yield was good. It's quite evident that the major products we obtained were circular 
DNA. The circular monomer products were then further purified and characterized with LC-MS and the 
expected molecular weights were observed (Figure 19b). 


2. RCA reaction: 

The ligation products were used for the RCA at different DNA concentrations. The products were 
characterized with agarose gel electrophoresis. From the gel, super long DNA single strand products 
were observed, which indicated the RCA reaction in solution was very successful. 
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Figure 20. 0.5% agarose gel electrophoresis analysis of the RCA products. 


3. RCA reactions on surfaces: 

From the solution RCA reaction results, it seems the products were much longer than lOkb. 
Therefore, if we use such long ssDNA to do surface modification, the surface coupling efficiency would 
be very low (active amino group is embedded). Besides, it would be difficult to get sufficient amount of 
long ssDNA products to do the surface modification. So, we decided to start with Approach B: to directly 
generate long ssDNA on surface. After surface RCA reaction, fluorescent probe DNA was added onto the 
surface. The probe DNA could only hybridize with the RCA product, but not the primer. The fluorescent 
image clearly showed that RCA reaction on the surface was successful. In addition, the negative sample 
gave very weak fluorescent signals, indicating a low background. 
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Figure 22. Fluorescent image of the surface RCA product incubated with fluorescent probe DNA. 


Currently, we are further studying the RCA reactions on PDMS surface and also testing the surface 
interaction mediated with the RCA products. 


Mid- to long-term challenges: 

Currently, we are focusing on the DNA-mediated mesoscale surface interaction. This is not a 
trivial problem, from both scientific and technical point of view. We envision there are a number of 
critical issues we are facing to: how to achieve flat, soft, and hydrophilic DNA encoded surface; how to 
minimize non-specific surface interaction; how to get good conformal contact between two surfaces; 
and how to optimize DNA hybridization between the surfaces. We are on the way to solving these 
problems. 

One of our mid-term challenges is that even if surfaces can interact through molecular 
recognition, it is quite another matter for them to adhere through transient contact while tumbling in 
solution. We might overcome this challenge by having a preliminary, weak nonspecific interaction 
between smartblock surfaces. This way, tumbling smartblocks would interact for a longer time, 
allowing the strong, specific interaction to occur, but allowing them to break up if mismatched. 

Furthermore, there are a number of kinetic issues which will need to be dealt with. Our scheme 
is on its face based on a thermodynamic argument: the correct arrangement of smartblocks is the 
lowest-energy one, because it maximizes the amount of surface interactions. Flowever, there are 
kinetic issues which may inhibit finding the thermodynamic minimum, such as off-center interactions 
and improper assembly order. 

The problem with off-center interactions is that, while the correct orientation maximizes the 
surface contact and is thus a thermodynamic minimum, it is more likely for tumbling smartblocks to 


26 





meet each other off-center. The best solution to this would be to lower the kinetic barrier for 
interacting smartblocks to "slide" against each other, allowing the thermodynamic minimum to be 
reached more readily. (Figure 23) 

Another kinetic issue is that if all smartblock surfaces are available to hybridize at once, the 
order of assembly will be uncontrolled, which could lead to problems like unoccupied lattice points 
(voids), and multiple nucleation sites causing regions of the object to partially coalesce separately. 
Preventing voids can be helped by selecting the best honeycomb lattice (e.g. truncated octahedra 
instead of cubes), but direct control of the assembly order would be a powerful solution to this problem. 
We have proposed a scheme for ordered self-assembly making use of DNA threshold gates which have 
been reported in the literature. (Figure 24) 


higher energy (less favorable) 



Figure 23. Schematic illustrating the problem of off-center interactions. The formation of the 
undesired off-center product is less thermodynamically favorable but more kinetically favored, 
indicating that smartblocks could get "stuck" in this arrangements. One way around this 
problem would be to minimize the kinetic barrier for directly converting off-center to on-center 
interactions, either through design or finding appropriate experimental conditions. 



Figure 24. Scheme for ordered self assembly. Left: Initially, only smartblocks in the "red" 
layer can assemble; "yellow" smartblocks in the next layer cannot assemble because their DNA 
handles are occluded by another oligonucleotide. The red smartblocks have handles which 
release a messenger oligonucleotide into solution upon assembly. Right: Once most of the red 
layer has assembled, the messenger oligonucleotide concentration becomes great enough that it 
triggers a threshold gate in solution that activates the yellow smartblocks by removing the 
occlusion on their handles. 
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Once these challenges have been overcome is complete, we may look into crosslinking 
chemistry or other secondary steps to supplement the DNA-mediated interactions in order to 
mechanically strengthen the object once it is assembled. (Milestones 3a-3b) 

We envision the following as possible downstream applications for this process (Phase III and 
beyond): 

• DNA aptamers that bind strongly to a small-molecule, protein, or nucleic-acid target can be readily 
evolved in the laboratory. The ability of aptamers to undergo hybridization changes upon target 
binding is well-established. This can be applied as programmable matter sensors that undergo 
desired structural changes in the presence of a target molecule (e.g. TNT, chemical or biological 
warfare agents, etc.) 

• The messenger DNAs released upon smartblock-megastaple hybridization, used in our ordered 
self-assembly scheme above, could also be used to induce the assembly or disassembly of other 
programmable matter structures by mediating or disrupting smartblock-megastaple pairing. In 
addition to allowing ordered assembly of one structure, it could enable communication between 
multiple structures (e.g. bistability), signal amplification, or communication with biological systems. 

• Our design for programmable matter may be amenable to directed evolution. Random libraries of 
megastaples induce the assembly of random three-dimensional structures. The structures are 
selected for a desired property (the ability to adopt a particular shape, mechanical strength under 
load, etc.). Structures that survive selection are disassembled and their megastaples are amplified, 
mutated, and reselected to evolve optimal programmable matter solutions to the task of interest. 


In conclusion, we believe that the properties of DNA make it very suitable as a platform for 
directing the assembly of programmable matter. We have achieved significant accomplishments 
towards this goal by developing the computational infrastructure needed to design DNA-mediated 
programmable matter, as well as creating and characterizing the components of our programmable 
matter scheme, the DNA-linked surfaces that will be used in smartblocks, the megastaples, and the 
super long repetitive single stranded DNA. We are currently testing these components to develop the 
first DNA hybridization-mediated mesoscale surface interaction, arguably the key tipping point of our 
proposal. 
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ABSTRACT : DNA nanotechnology provides a versatile foun¬ 
dation for the chemical assembly of nanostructures. Plasmonic 
nanoparticle assemblies are of particular interest because they 
can be tailored to exhibit a broad range of electromagnetic 
phenomena. In this Letter, we report the assembly of DNA- 
functionalized nanoparticles into heteropentamer clusters, 
which consist of a smaller gold sphere surrounded by a ring of 
four larger spheres. Magnetic and Fano-like resonances are 
observed in individual clusters. The DNA plays a dual role: it 
selectively assembles the clusters in solution and functions as an insulating spacer between the conductive nanoparticles. These 
particle assemblies can be generalized to a new class of DNA-enabled plasmonic heterostructures that comprise various active and 
passive materials and other forms of DNA scaffolding. 

KEYWORDS: Plasmonics, DNA, colloidal self-assembly, nanoshell, magnetic dipole, Fano resonance 



S ub-wavelength-scale metallic particles are a basis for nano¬ 
scale light manipulation because they support localized sur¬ 
face plasmon resonances, which are oscillations of free electrons 
in metal that couple with electromagnetic waves. 1 By synthesiz¬ 
ing particles into specific shapes and engineering their assembly, 
it is possible to construct nanoscale chemical sensors, 2 plasmonic 
rulers, 3 optical nanocircuits, 4 and metamaterials. 5,6 DNA nano¬ 
technology 7 is a vehicle for the controllable assembly of nano¬ 
particles 8 because it enables the positioning of particles with 
nanoscale precision and the tailoring of their binding interac¬ 
tions. While simpler implementations of DNA particle assembly 
involve controlled nanoparticle aggregation, 9,10 other efforts 
have focused on the construction of well-defined clusters and 
lattices. For example, micrometer-scale dielectric particles have 
been assembled into tetrahedral and octahedral clusters, 11 and in 
other schemes, trimer clusters, 12,13 tetrahedral clusters, 14 chiral 
helical assemblies, 15 and two- and three-dimensional lattices 16,17 
of plasmonic nanoparticles have been constructed. It is tantaliz¬ 
ing to envision these nanostructures as useful optical structures; 
however, these plasmonic assemblies have yielded little optical 
data in the literature because they typically utilize very small 
metallic particles (diameter <20 nm). With such small particles, 
optical measurements on individual nanostructures become 


extremely difficult due to their small scattering cross sections. 
Larger optical signals can be obtained from ensemble measure¬ 
ments, but these suffer from sample heterogeneities that can 
weaken or completely eliminate the observation of certain optical 
resonances. 18 Another issue with small particles is that, unlike 
large particles, they do not support many plasmonic modes; 
higher order resonances require retardation for their excitation, 
and Fano-like resonances, which will be discussed later, require 
strong linewidth broadening from radiative damping and retar¬ 
dation that exist only in large particle systems. 

In this study, we use DNA to construct heteropentamer 
clusters consisting of a solid gold nanosphere surrounded by 
four larger nanoshells. 19 The DNA route to particle assembly 
has distinct advantages over previous random capillary assembly 
methods: 6 it introduces more specific and programmable inter¬ 
actions between metallic particles of different sizes and types and 
provides a potential route to more sophisticated three-dimen¬ 
sional particle assembly. The construction of the heteropentamer 
was partially motivated by prior calculations demonstrating 
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Figure 1. DNA-mediated assembly of plasmonic heteropentamers. (a) To assemble the clusters, gold nanospheres and nanoshells are functionalized 
separately with thiolated DNA strands, which form self-assembled monolayers on the particle surfaces. They are then mixed and incubated together at 
room temperature, where they assemble into "loosely” bound pentamers. Finally, the clusters are dried on a hydrophilic substrate, where capillary forces 
compress the pentamers into a close-packed two-dimensional configuration, (b) The nanospheres are functionalized with partially dsDNA molecules 
consisting of three regions: a gold attachment region comprising a thiol group and 5-mer polyT sequence, a 75bp ds "spacer” region, and a 20-mer 
"linker” region. Prior to nanosphere attachment, the thiolated 100-mer DNA (red) is mixed with complementary 75-mer strands (green) to form rigid 
dsDNA segments that effectively increase the nanosphere size in solution (dotted circle). The nanoshells are functionalized with ssDNA (blue) 
comprising a 50-mer polyT spacer sequence, followed by a 20-mer linker that complements the nanosphere linker. 


Fano-like resonances for this cluster geometry; 20 sufficiently large 
nanoparticles are used here so that the Fano minima are clearly 
observable. The ratio of the nanoshell diameter to nanosphere 
diameter is set to approximately 2.4 to ensure a close packed 
geometry and strong coupling between all neighboring particles. 

The self-assembly process is outlined in Figure la. First, the 
nanospheres and nanoshells are functionalized separately with 
different thiolated DNA molecules, 21 which form self-assembled 
monolayers on the particle surfaces. The outermost 20-mers of 
DNA on the nanospheres and nanoshells are complementary, 
thus facilitating specific nanosphere—nanoshell binding while 
minimizing interactions between particles of the same type 
(Figure lb). The particles are then cleaned and incubated to¬ 
gether at high salt concentration and room temperature, where 
loosely packed pentamers are formed in three dimensions in 
aqueous suspension. The ratio of the number of nanoshells to 
nanospheres is set to 12:1 to enhance pentamer yield while 
limiting the assembly of large nanoshell—nanosphere aggregates, 
which were observed for particle mixtures with smaller ratios. 
Finally, the clusters are air-dried on a hydrophilic substrate at 
room temperature. During the drying process, capillary forces 
compress these loosely packed three-dimensional structures 
into two-dimensional pentamers. Both single-stranded DNA 
(ssDNA) and double-stranded DNA (dsDNA) collapse and 
pack to yield dense 2 nm thick dielectric spacers between all 
neighboring nanoparticles, ensuring their strong and controlled 
optical coupling. As a control experiment, nanoparticles were 
functionalized with non-complementary DNA sequences and 
mixed together, and no pentamers were observed. It is difficult to 
quantify the overall yield of these pentamer clusters (further 
discussed in the Supporting Information); nevertheless, it is easy 
to find many pentamers on a TEM grid, which is sufficient for 
single nanostructure experimentation. We also found other 
DNA-enabled clusters, ranging from dimers and trimers to small 
aggregates. Heteropentamer yields may be further enhanced by 


increasing the number ratio of nanoshells to nanospheres and 
employing separation methods, such as density gradient centri¬ 
fugation, to distill pentamers from single nanoshells and other 
cluster types, and this will be the subject of future study. 

The assembly of pentamers in solution involves the sequential 
attachment of nanoshells onto a nanosphere (Figure 2a). The 
first association between a single nanoshell and nanosphere is 
straightforward: Brownian motion brings the two particles into 
close proximity and their surface-attached DNA hybridizes. 
However, the association of additional nanoshells onto a nano¬ 
sphere becomes more difficult kinetically for two principal 
reasons (Figure 2b). The first involves steric hindrance: already 
associated nanoshells will physically block other nanoshells from 
getting close to and associating with the nanosphere. The second 
is due to the lack of hybridizable DNA on the nanosphere surface. 
Single-stranded DNA is a polymer with ~1—2 nm persistence 
length in solutions containing high salt concentration (ref 22), 
such that multiple DNA strands from a single nanoshell can 
attach to multiple nanosphere DNA strands in a polyvalent inter¬ 
action. In our system where the length of DNA is on the order of 
the nanosphere diameter, the DNA from two or three nanoshells 
can associate with most of the nanosphere strands, leaving very 
few free nanosphere strands for additional association events. In 
order to overcome these problems, we design the 75 base double- 
stranded spacer on the nanosphere. Since the persistence length 
of dsDNA is 50 nm, this rigid spacer effectively increases the size 
of the nanosphere in solution by 50 nm 23 (Figure lb), alleviating 
the steric hindrance problem. Also, the rigidity of the dsDNA and 
relatively short length of the ssDNA linker reduce the polyvalent 
association between the nanoshells and nanospheres. 

Individual heteropentamers are identified using transmission 
electron microscopy (TEM), and scattering spectra from indivi¬ 
dual clusters are measured in the near-infrared frequency range 
using dark-field microscopy. (See Supporting Information for 
setup.) Here, the incident light is s-polarized (i.e., the electric 
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Figure 2. Detailed view of plasmonic cluster formation in solution, (a) The pentamer assembles in solution by the sequential attachment of nanoshells 
onto a nanosphere. Successful pentamer assembly requires each of these intermediate steps to be stable (i.e., the nanoshells remain attached) and 
kinetically accommodating to additional nanoshell association events, (b) There exist at least two issues that can impede the attachment of additional 
nanoshells to a cluster. One is that already associated nanoshells can sterically hinder other single nanoshells from getting close to and attaching to the 
nanosphere (left). This is especially a problem if the nanosphere DNA strands are too short. Another is that if the nanosphere DNA is too long and 
floppy, DNA from just two or three nanoshells can associate with most of the nanosphere DNA, reducing the total number of available linker DNA on the 
nanosphere and suppressing additional nanoshell association events (right). Both of these issues are addressed by utilizing the long dsDNA spacer on the 
nanosphere. 





Figure 3. TEM images and scattering spectra of individual heteropentamer and quadrumer clusters, (a) s-polarized spectra of a pentamer for different 
electric field polarization angles. Both the experimental and theoretical spectra are characterized by a broad electric dipole peak and a narrow Fano dip 
near 1000 nm. The in-plane isotropy of these spectra is consistent with the symmetry of these clusters. The simulated geometry is based on the TEM 
image and uses a nanosphere diameter d = 74 nm, nanoshell [r 1; r 2 ] = [62.5, 92.5] nm, and interparticle gaps of 2 nm that are filled with dielectric, (b) 
TEM image and s-polarized spectra of a symmetric quadrumer. These spectra exhibit a broad electric dipole peak but no clear Fano dip, indicating that 
the nanosphere is necessary for this cluster system to exhibit a strong Fano-like resonance. The experimental spectra generally match those simulated 
using the same nanoshell geometry and gap dielectric as in (a). 


field is in the plane of the cluster) and the polarization angle of 
the electric field relative to the nanostructure is controlled. 
Scattered light from an individual structure is selected with a 
confocal-type technique, by placing a pinhole at a reimaged focal 
plane and allowing light from only the structure of interest to 
enter the spectrometer. Spectra of a single pentamer for three 
different polarization angles are shown in Figure 3a and are 


characterized by a broad electric dipole resonance spanning the 
entire range of the plot. There exists a narrow and asymmetric dip 
near 1000 nm, which is a Fano resonance. These resonances arise 
from interference between a superradiant “bright” mode and 
a subradiant “dark” mode and have been studied extensively 
in plasmonic systems. 5 ' 24 These experimental spectra match 
those calculated using the finite element simulation software 
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Figure 4. Theoretical analysis of the bright and dark modes involved with Fano-like resonances, (a) Extinction spectrum and surface charge plots of the 
heteropentamer excited at normal incidence with a 45° polarization angle. The charge density of the bright mode, which is peaked at 1050 nm (pink 
dashed line), shows the charge oscillations on each nanoparticle oriented in the same direction, yielding a large cluster dipole moment. The dark mode at 
980 nm (black dashed line) shows the charge oscillations on each nanoparticle oriented in different directions, yielding a small cluster dipole moment 
and suppressed radiative losses. There is clear capacitive coupling between the nanosphere and the two nanoshells above and below the nanosphere, 
which redshifts the mode close to the bright mode peak, (b) Extinction spectrum and surface charge plots of the quadrumer excited with the same 
conditions as in (a). The bright mode peak is near 1040 nm (pink dashed line) and the charge distribution shows a large cluster dipole moment. The dark 
mode peaked at 650 nm (black dashed line) shows the total cluster dipole moment is small. The absence of capactive coupling between the nanoparticles 
prevents strong overlap in frequency between the dark and bright modes. 


COMSOL (Figure 2a) and are nearly independent of incident 
polarization angle due to the D 4 ^ group symmetry of the cluster, 
which supports isotropic in-plane resonances in the quasistatic 
limit. 

To further probe the interaction of the nanosphere with the 
nanoshells in the pentamer, we identified and optically charac¬ 
terized symmetric nanoshell quadrumers. These nanostructures 
assembled by entirely different means from the pentamer: here, 
four single nanoshells trapped in a drying droplet clustered 
together by capillary forces, in similar fashion to clusters as¬ 
sembled previously. The strong electrostatic repulsion between 
the nanoshells, provided by the negatively charged DNA, helps 
them arrange in a non-close packed configuration; it is noted that 
in previous experiments utilizing uncharged PEG-functionalized 
nanoshells, symmetric quadrumers were not found. 6 The spectra 
from an individual quadrumer display a broad electric dipole 
resonance but no strong Fano dips, indicating that the nanosphere 
in the pentamer is required for strong Fano-like resonances to be 
observed (Figure 3b). As with the pentamer, this cluster possesses 
isotropic in-plane resonances due to its group symmetry. 

The Fano-like resonances in these clusters are further analyzed 
by examining the surface charge distributions of their bright and 
dark modes (Figure 4). In both the quadrumer and pentamer, the 
bright modes are characterized by nanoparticle polarizations 
oriented in the same direction. As such, the cluster dipole mo¬ 
ments are large and the mode strongly redshifts due to strong 
capacitive coupling between neighboring nanoparticles. Physi¬ 
cally, this coupling arises from the attractive quasi-static interac¬ 
tion between surface charges on these nearly touching particles. 25 
The dark modes of each cluster, however, have very different 
spectral positions. The pentamer dark mode charge distribution 
shows that the small nano sphere capacitively couples with two of 
the nanoshells, which strongly redshifts the mode. The total 
dipole moment of this cluster is small but non-negligible, which 
limits the magnitude of the Fano dip. The quadrumer dark mode, 
which has been studied elsewhere, 26 exhibits little capacitive 
coupling between adjacent nanoshells. As a result, the mode is 



Wavelength (nm) Wavelength (nm) 



Figure 5. Magnetic dipole resonances in nanoparticle clusters, (a) 
Experimental and theoretical cross-polarized spectra of the pen¬ 
tamer at 0° polarization angle (left) reveal a narrow peak near 
1400 nm, which matches the peak position of the calculated out-of- 
plane magnetic dipole moment of the cluster (right). The inset is a 
quasistatic mode plot of the surface charges (color) and displace¬ 
ment current (arrows) of the magnetic dipole mode, (b) Experi¬ 
mental and theoretical cross-polarized spectra of the quadrumer 
with the same electric field polarization angle (left) also reveal a 
narrow magnetic dipole peak near 1400 nm that matches the peak 
position of the calculated magnetic dipole moment (right) of the 
cluster. The inset is a quasi-static mode plot of the magnetic dipole 
mode and shows a circulating current around the ring of nanoshells 
similar to that in (a). These plots and calculated magnetic dipole 
moments show that the magnetic dipole mode in the pentamer and 
quadrumer is due to near-field interactions between the nanoshells 
and that the mode in the pentamer is effectively decoupled from the 
nanosphere. 
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blueshifted beyond the bright mode and there is no clearly visible 
Fano minimum. As observed with other nanocluster structures 
previously studied/ the key to engineering a strong Fano dip is to 
design bright and dark modes that exhibit similar levels of 
capacitive coupling between particles, thereby enforcing strong 
spectral overlap between the two modes. 

We also analyze magnetic dipole modes in these clusters. 
These resonances are excited by the magnetic component of the 
incident electromagnetic field and were previously measured in 
nanoshell trimers. 6 They are predicted to be supported here by 
the outer ring of nanoshells, which form a closed loop of metallic 
nanoinductors and dielectric nano capacitors 27 in similar fashion 
to the trimer. These modes are not clearly visible in the scattering 
spectra in Figure 5 because they weakly scatter compared to 
the electric dipole; a cross-polarizer oriented 90° relative to 
the incident light polarization can be placed after the collec¬ 
tion objective to filter out elastically scattered electric dipole 
radiation. 6 The cross-polarized spectra of the quadrumer and 
pentamer are shown in Figure 5 and they both exhibit clear, 
narrow peaks near 1400 nm. The positions of these peaks match 
those of the calculated magnetic dipole moments, confirming 
that they are magnetic resonances. There is less background in 
these spectra compared to those of the trimer previously studied 
because unlike the trimer, these clusters have inversion symmetry 
and, based on the criteria outlined in ref 28, do not support 
extensive optical activity. The simulated electric field profile and 
displacement current of this mode in both cluster types show a 
circulating current around the quadrumer ring, which is a hall¬ 
mark of a classical magnetic dipole. The matching magnetic 
dipole moments from both clusters indicate that the presence of 
the nanosphere in the pentamer does not affect the magnetic 
resonance from the nanoshell ring. This can be further under¬ 
stood by group theory, in which the irreducible representation of 
the pentamer can be expressed as T Pent = r Nanosphere + r auad . 
The magnetic mode in the quadramer has an irreducible 
representation of A 2g while the nanosphere has a irreducible 
representation of E lu ; as modes with different representations do 
not couple with each other, the nanosphere does not affect the 
magnetic mode. 

The DNA-enabled assembly of nanoparticles can be general¬ 
ized to a broad range of two- and three-dimensional hetero¬ 
clusters. These are not limited to passive plasmonic particles and 
can include other types of dielectric, nonlinear, active, and 
organic materials to create new functional nanostructures such 
as active antennas, 29 surface plasmon lasers, 30 clusters with 
tailorable hot spots, 31 and metamaterial fluids. 32 Other forms 
of DNA nanotechnology, such as DNA origami, 33 have great 
potential as rigid scaffolds with high spatial resolution for particle 
assembly. In this study, the pentamers were assembled in two 
steps, first by DNA to “loosely” assemble the clusters in suspen¬ 
sion and then by capillary forces to “compact” the clusters in two 
dimensions; loosely assembled clusters can be compacted or 
manipulated in three dimensions in suspension via depletion 34 or 
optical forces 35 or by tailoring new chemical interactions be¬ 
tween particles, 36 and new regimes of capillary assembly can be 
employed using patterned substrates, 37 surfactants, or other 
additives to further control cluster assembly on substrates. The 
merging of biomaterials like DNA with plasmonic nanostructures 
is suggestive of new forms of plasmon-enhanced biomolecular 
detection schemes, dynamically reconfigurable nanostructure 
geometries, and even direct integration and assembly of na¬ 
noclusters within biological systems. 
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ABSTRACT : We developed a sequential strand-displacement 
strategy for multistep DNA-templated synthesis (DTS) and 
used it to mediate an efficient six-step DTS that proceeded in 
35% overall yield (83% average yield per step). The efficiency of 
this approach and the fact that the final product remains linked 
to a DNA sequence that fully encodes its reaction history 
suggests its utility for the translation of DNA sequences into 
high-complexity synthetic libraries suitable for in vitro selection. 


D NA-templated synthesis (DTS), the use of DNA hybridiza¬ 
tion to dramatically increase the effective molarity of reactants 
linked to oligonucleotides, is a powerful strategy to control chemical 
reactivity in a DNA sequence-programmed manner. 1-4 Because the 
reaction products of DTS are encoded by the sequences of the 
associated DNA templates, they can be subjected to in vitro selection 
followed by PCR amplification and DNA sequence analysis to enable 
the discovery of functional small molecules, 3,5-7 synthetic poly¬ 
mers, 8- 10 or novel chemical reactions. 11-15 We recendy reported the 
three-step DNA-templated synthesis of a 13,824-membered small- 
molecule macrocycle library. 16 The library was subjected to in vitro 
selection for binding affinity to a variety of proteins of biomedical 
interest, ultimately yielding a new class of macrocyclic kinase inhi¬ 
bitors. 6 Other complementary approaches to generating DNA- 
encoded libraries have led to the discovery of bioactive small 
molecules, 17-29 including a number of examples in the pharmaceu¬ 
tical industry. 5,30,31 

Generating DNA-encoded small molecules of significant 
structural complexity requires multistep DNA-programmed or 
DNA-tagged synthesis. 3,5, 3,30,32-37 A number of strategies have 
been developed to enable multistep DTS. The simplest uses a 
DNA template strand containing several codons and reagents linked 
to complementary anticodon oligonucleotides that are added succes¬ 
sively. 3,32-34 While this approach is conceptually straightforward, it 
requires several manipulations after each step that increase required 
time and effort and can decrease overall yields. The relative geometry 
between reactants on the template and reagent strands also 
changes after each step in this approach, potentially altering reaction 
efficiencies. 38 More complex self-assembled DNA structures and 
devices can also mediate multistep DTS. For example, a DNA three- 
way junction that contains multiple reagents at the junction has been 
developed for the construction of DNA-encoded peptides. 5 We 
developed a DNA mechanical device that moves along a DNA track 
and mediates autonomous multistep organic synthesis in a single 
isothermal solution. 35 McKee et al. recendy used a DNA strand- 
exchange strategy to achieve a three-step DTS in which products are 


swapped between new and old DNA strands with the assistance of a 
“remover strand” that displaces expended reagent oligonucleotides. 36 

Despite these significant advances and the diversity of ap¬ 
proaches to generating multistep DTS products, all multistep 
DNA-templated small-molecules syntheses reported to date have 
used only three or fewer DNA-templated steps, and overall yields 
are generally low (typically <10%). Here we present a new 
strand-displacement strategy for multistep DTS and its use to 
mediate a six-step synthesis with an overall yield of 35% (average 
yield of 83%). By providing products of six-step DNA-pro¬ 
grammed reaction sequences in good overall yield, the approach 
presented here may provide access to high-complexity DNA- 
templated small-molecule libraries. 

Our strategy exploits “toehold displacement,” the known 
ability of a single-stranded DNA oligonucleotide (AB) to invade 
an asymmetric DNA duplex (A'B'iB) that contains a single-stranded 
hybridization site (A r ) for the invading strand. 39 Once A:A! hybri¬ 
dization takes place, base pairing with the invading strand continues, 
ultimately resulting in strand displacement of the shorter, and 
therefore less-favorably hybridized, B strand. Displacement results 
in the formation of a new Watson—Crick complex suitable for DTS. 
We hypothesized that this approach could represent a highly efficient 
and very simple way to access products of several consecutive DNA- 
programmed reactions while preserving the correspondence between 
DNA sequence and reaction product structure that is required for in 
vitro selection. 

The application of this “toehold displacement” strategy to 
multistep DTS is summarized in Figure 1. A single-stranded 
DNA template (T) consists of a 16-base initial reaction site 
(black) followed by five consecutive 8-base coding segments 
(colored) that also serve as toeholds to initiate sequential DNA 
strand displacement. In the first step, substrate DNA SI, which is 
tethered to the first reactant, hybridizes to T, initiating DTS. If 
self-cleaving reagents are used, the first reactant group is trans¬ 
ferred from SI to T as a natural consequence of the DNA- 
templated reaction. After sufficient time to react has elapsed, the 
second substrate strand S2 is added without requiring product 
purification or any other manipulation. S2 partially hybridizes 
with T through the red toehold region and displaces SI by 
forming a more stable DNA duplex due to the fact that the 
oligonucleotide in S2 is eight nucleotides longer than that in SI. 
As a result, SI spontaneously leaves the Watson—Crick complex 
and the second reactant's group is transferred to the T-linked 
intermediate through a DNA-templated reaction. Thereafter, an 
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Figure 1. Toehold displacement—mediated multistep DNA-templated 
synthesis. In each step, hybridization between the template and substrate 
oligonucleotides triggers a DNA-templated reaction. The template 
strand contains multiple toehold regions (colored) that allow successive 
docking of new substrates and DNA strand displacement of oligonu¬ 
cleotides following each reaction without requiring purification after 
each step. 


even longer substrate strand S3 is added to displace S2 and 
initiate a third reaction. 

So long as unpaired toeholds (colored regions in Figure l) 
remain in T, this process can be iterated to generate successively 
longer or more complex products without requiring any intervention 
during the reaction sequence beyond the simple addition of 
successive substrates to the same reaction vessel. Importantly, the 
effective molarity of the reactants during successive DNA-templated 
steps in this strategy remains relatively constant since the reactive 
groups for each step remain at the end of the double-stranded 
region. 38 As a result, we anticipated that reaction yields would remain 
higher over the course of many steps than with traditional multistep 
DTS methods. Crucially, the DNA templates in this strategy retain all 
product-encoding information throughout the DTS process, thus 
enabling the multistep products to be subjected to in vitro selections 
for desired properties, in contrast with a different approach to 
multistep synthesis that uses DNA strand exchange. 36 

We designed and executed a six-step DNA-templated reaction 
sequence using this approach. The DNA template (T) contains a 
3' primary amine group, and the first five substrates (SI—S5) 
contain S' N-hydroxysuccinimidyl (NHS) ester-linked amines. 32,35 ' 37 
The oligonucleotide—reagent bond is automatically cleaved as a 
result of the amine acylation reaction, thus enabling spontaneous 
transfer of the amino acid building block from each substrate to the 
template. The transferred amino acid displays another primary amine 
group that is poised to further react with the next NHS ester to 
continue the multistep DTS. The last substrate (S6) contains a biotin 
NHS ester and lacks a primary amine group (Figure 2). 

In a typical experiment, T (l juM) and SI (1.2 //M) were 
combined in aqueous buffer containing 50 mM MOPS and 10 mM 
Mg(OAc) 2 at pH 7.5. After 30 min at room temperature, the second 
substrate S2 was added into the system, followed by another 30-min 



Figure 2. (a) Reactant structures; (b) expected six-step DNA-tem¬ 
plated reaction product. 

incubation at room temperature. Substrates S3—S6 were similarly 
added at successive 30-min intervals. During the entire process, the 
temperature and buffer conditions were constant. The six-step reac¬ 
tion was therefore carried out in a single vessel at room temperature 
over 3 h. The crude reaction mixture was desalted by gel filtration and 
analyzed by high-resolution ESI mass spectrometry. For comparison, 
reactions containing only T, T+Sl, T+S1+S2, T+S1+S2+S3, 
T+S1+S2+S3+S4, T+S1 + S2+S3+S4+S5, and Sl-hT,^ (a 
template containing a scrambled DNA sequence) were also sub¬ 
jected to the reaction conditions and analyzed by high-resolution 
mass spectrometry. 

The unreacted template and all six expected reaction products 
corresponding to one-, two-, three-, four-, five-, and six-step reaction 
sequences were observed with high mass accuracy (Figure 3). Each 
reaction resulted in highly efficient formation of the desired multistep 
product, such that the predominant template-linked ion observed 
from each reaction corresponded to the expected one- to six-step 
reaction product (Figure 3). The observation that in the six-step 
reaction mixture the template-linked full-length hexaamide product 
provided the strongest m/z signal among all template-linked ions 
suggested that the six-step product might be formed in excellent 
overall yield. The control reaction containing the mismatched 
template and SI did not generate detectable product (Figure 3), 
indicating that hybridization with the template is required to effect 
product formation under these conditions. 

To quantitate reaction yields during all six reaction sequences, 
we used denaturing polyacrylamide gel electrophoresis (PAGE). 
Because each building block adds significant mass to the template, 
the full-length product and all the intermediate products migrate 
differendy during PAGE. To achieve better separation during PAGE, 
all reaction products were digested with restriction endonuclease 
Dde I, which cleaves the 56-nt template into a 22-nt 3' fragment 
containing the reaction products. To enable accurate quantitation of 
reaction yields, a single fluorescein group was incorporated into the 
template (prior to DTS) near the 3' end, allowing the 22-nt frag¬ 
ments to fluoresce with equal intensity. Densitometry of fluorescein- 
visualized PAGE gel (Figure 4) revealed that the three-step reaction 
sequence (lane 4) generated the corresponding three-step product in 
55% overall yield, and the six-step reaction sequence (lane 7) 
produced full-length hexaamide product in 35% overall yield, corre¬ 
sponding to ~83% average yield per step (see the Supporting 
Information [SI] for additional details). An independent experiment 
to determine reaction yields using nonfluoroscein-labeled T in which 
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Figure 3. Mass spectrometric analysis of multistep DNA-templated 
reaction products generated by sequential strand displacement. Eight 
reactions using the species shown were separately performed as de¬ 
scribed in the text. The mass range shown includes all template-linked 
species. Observed m/z values for ions of —5 charge are shown in black; 
expected m/z values are shown in red. The arrow in the last spectrum 
indicates the expected m/z value of the product corresponding to 
reaction between SI and a scrambled template (T mis ). 


Ddel-digested reaction products were analyzed by PAGE followed by 
staining resulted in 38% overall yield of six-step reaction product (see 
the Si). To our knowledge, these results represent both the highest 
yielding and longest multistep DNA-templated small-molecule 
syntheses to date. 32,35-37 

Applying this strategy for future multistep DNA-templated 
library syntheses will require that reagent-linked oligonucleotides 
each be prepared as mixtures containing all possible anticodons 
for prior DNA-encoded steps. The recent availability of inexpen¬ 
sive custom oligonucleotide pools from chip-based parallel 
synthesis 40-43 provides practical access to these materials. In 
addition, the fact that the specificity-determining step of each 
reaction is hybridization of the 8-base toehold, which encodes the 
current reagent, implies that mismatch hybridization elsewhere 
in the template:reagent oligonucleotide duplex will not affect the 
fidelity of the multistep DTS process. Moreover, longer DNA 
templates could conceivably give rise to products of more than six 
DNA-programmed reactions. Indeed, we have observed that a 
64-base, 80-base, and 96-base DNA strand in a duplex, sufficient 
to encode up to ~10 DNA-templated reactions, can be effec¬ 
tively displaced using a 72-mer, 88-mer, or 104-mer (see the SI 
and Figure S3). 
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Figure 4. Denaturing PAGE analysis of the multistep DTS reactions 
analyzed in Figure 3. Each template-linked species contained a single 
fluorescein group to enable uniform visualization. The reaction products 
were digested with the endonuclease Ddel prior to PAGE analysis to 
cleave the product-linked oligonucleotides to 22-mers. Lanes 1 and 8: 
unreacted template; lanes 2—7: DNA template after successive reaction 
with one, two, three, four, five, or six reagents, respectively. 


In summary, we developed a rapid, one-pot, highly efficient 
multistep DNA-templated synthesis strategy that is driven by 
successive toehold DNA strand displacement. The strategy was 
used to generate a six-step DTS reaction product in 35% overall 
yield (83% average yield per step), such that the desired product 
was the major observed template-linked species after six steps 
with no purification or reaction manipulation required beyond 
addition of each substrate every 30 min. Since the reaction 
history of the products remains encoded by their associate 
templates, this strategy may enable future efforts to generate 
high-complexity (e.g., n 6 ) DNA-templated libraries suitable for in 
vitro selection. 
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